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Abstract 

We review recent numerical work investigating the equilibrium phase diagram, and the dynamics, 
of the cholesteric blue phases. In equilibrium numerical results confirm the predictions of the classic 
analytical theories, and extend them to incorporate different values of the elastic constants, or the 
effects of an applied electric field. There is a striking increase in the stability of blue phase I in 
systems where the cholesteric undergoes helical sense inversion, and the anomalous electrostriction 
observed in this phase is reproduced. Solving the equations of motion allows us to present results 
for the phase transition kinetics of blue phase I under dielectric or flexoelectric coupling to an 
applied electric field. We also present simulations of the blue phases in a flow field, showing how 
the disclination network acts to oppose the flow. The results are based on the Landau-de Gennes 
exapnsion of the liquid crystal free energy: that such a simple and elegant theory can predict such 
complex and subtle physical behaviour is remarkable. 

PACS numbers: 61.30.Mp, 64.70.mf, 83.80.Xz, 61.30.Gd 



I. INTRODUCTION 



One of the spectacular successes of de Gennes' Landau theory of hquid crystals is that 
it captures, without need for extension or modification, many of the remarkable and subtle 
phase transitions and properties of liquid crystals. Blue phases provide an especially vivid 
example of this. The general features and properties of these unique mesophases are both 
'beautiful and mysterious' [Ij and yet they can be understood qualitatively and quantita- 
tively using de Gennes' Q-tensor theory of liquid crystal phase transitions. 

Blue phases are found in highly chiral liquid crystals between the high temperature 
isotropic fluid and the lower temperature cholesteric phases [U El El IH El El E] • They are 
remarkable mesophases, exhibiting a brightly coloured texture of individual, micron sized 
platelets. The bright colour indicates selective reflection due to a periodic structure, much 
like in an ordinary crystal, but with a much larger characteristic length scale, and indeed 
the reflection spectra show Bragg peaks that can be indexed by cubic space groups with 
lattice constants of several hundred nanometers. Furthermore, individual platelets of mon- 
odomain crystals themselves show distinctive facetting corresponding to Miller planes of the 
lattice structure and with the faces growing in a sequence of steps |8]. Yet, blue phases are 
not crystals in the traditional sense: they have no long range positional order and are full 
three-dimensional fluids. The crystalline order is in the orientational degrees of freedom of 
the liquid crystal. 

The key to understanding the properties of the blue phases was in realising that the 
locally preferred order in the chiral liquid crystal is one of double twist, with the molecules 
adopting helical ordering along two, perpendicular opposed to the usual single twist 

of an ordinary cholesteric helix [9l [TOl [HI [12]. However, it is not possible to construct a 
global state with helical ordering in two directions, without introducing disclination lines 
into the structure. Therefore the blue phases consist of local, cylindrical, regions of double 
twist separated by a regular disclination network. 

A tutorial, two-dimensional example of a blue phase texture is shown in Fig. [l} This 
demonstrates how local regions of double twist can be pieced together only by introducing a 
square array of topological defects. In physical, three-dimensional, crystals the disclination 
structures are more complex; three distinct blue phase textures have been identified upon 
cooling from the isotropic liquid. Typical experimental phase diagrams are shown in Fig. |2} 
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Blue phases I and II exhibit textures with cubic symmetry corresponding to the space groups 
O^^ and respectively [7j. These are illustrated in Fig. [s] which shows both the disclination 
networks and the topology of the lattice of double twist cylinders in the two phases. Blue 
phase III is less well understood; it has an amorphous structure with the same symmetry as 
the isotropic fluid ^E^. 

The unique combination of crystalline order with lattice constants comparable to the 
wavelength of visible light and full three dimensional fluidity make the blue phases ideal 
for technological uses such as fast light modulators, photonic crystals or tunable lasers [13, 
\W\ [T71 US]- The principal obstruction to their application was the very limited temperature 
range, ~ IK, over which the phases were stable. However, recently, the stability range 
has been extended to as much as 60K, including room temperature, by the addition of 
bimesogenic molecules or photo-crosslinking of polymers [191 ED] • This has opened the way 
to device applications and, in May 2008, Samsung Electronics unveiled a prototype 'blue 
phase mode LCD' [211 [22] • 

The theoretical description of the general features and properties of the blue phases was 
developed in the 1980s in a series of important papers, summarised in the review [7]. Two 
approaches were followed, a low chirality theory, based on the Frank director field descrip- 
tion [9l [T0| and a high chirality theory using the Landau-de Gennes Q-tensor [23l [2l] . The 
theoretical approaches were able to correctly predict the symmetries of the blue phases and 
to give a good account of the phase diagram, although with some discrepancies. For exam- 
ple, additional textures with space groups O^^ and were found to be stable within the 
parameter range corresponding to the experiments, but were not observed experimentally. 

The renewed experimental interest in the blue phases provided the motivation to revisit 
them theoretically. Since the original analytic investigations, computers have become a 
powerful tool for studying complex fluids, and it is now feasible not only to revisit and 
extend the calculations of the blue phase phase diagrams [231201, but also to obtain results 
on blue phase kinetics and hydrodynamics [271 [28] . The aim of this article is to review 
recent research showing how numerical simulations of the blue phases are allowing us to 
gain further insights into their behaviour. 

In Section II we describe the Landau-de Gennes equations that have proved so successful 
in modelling the blue phases and outline how they can be solved numerically. In particular 
we discuss how the size of the unit cell must be allowed to vary to correctly minimise the 
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free energy. Then, in Section III, we revisit the classic calculations of the thermodynamic 
phase diagram showing the changes that result when the approximations necessary to make 
the analytic calculations feasible can be circumvented by a numerical solution. We also 
show that including additional terms in the Landau-de Gennes free energy can significantly 
change the phase diagram. Most dramatic among the results is the increase in the stability 
of blue phase I obtained in systems where the cholesteric undergoes helical sense inversion. 

Many of the potential technological applications of the blue phases rely on the response 
of the material to an electric field. In the blue phases the effects of an electric field are 
particularly rich, including continuous distortions in the size and shape of the unit cell [291 
[301 [31] and a series of field-induced transitions to new blue phase structures, not stable in 
zero field [32], [33l [Ml [35] . In Section IIIC we show how the Landau-de Gennes expansion 
correctly predicts electrostriction and, in particular, the anomalous electrostriction of blue 
phase I. 

We next, in Section IV, summarise the equations of motion which describe blue phase 
hydrodynamics. These can be used to investigate the kinetics of transitions between blue 
phases, and their response to an applied electric field or to an imposed flow. Section V is 
devoted to examples demonstrating the application of the dynamical equations. We describe 
an investigation into the viscoelastic properties of the blue phases as they are subject to a 
Poiseuille flow field. We then describe the kinetics of the phase transitions induced in blue 
phase I as it is coupled to an electric field, either via a dielectric, or a fiexoelectric, term in 
the free energy. This allows us to propose a candidate structure for blue phase X. 

II. MINIMISING THE LANDAU-DE GENNES FREE ENERGY 
A. Landau expansion in terms of the Q-tensor 

One of de Gennes' vital contributions to the theory of liquid crystals was to identify the 
Q tensor as a suitable order parameter pLj . A Landau expansion of the free energy of liquid 
crystals in terms of Q has proved useful as a starting point for interpreting liquid crystal 
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phase diagrams [B [3 |23l [261 136] 




Here Aq is a constant with the dimensions of an energy density, 7 plays the role of an 
effective temperature for thermotropic liquid crystals, go defines the helical pitch and the 
Lij are elastic constants. The expansion, taken to second order in the derivatives of Q, can 
only account for two independent Frank elastic constants and imposes that the magnitude 
of splay and bend are equal [7]. In order to remove this constraint, and to allow for a 
temperature dependent helical pitch, it is necessary to consider, at least, terms cubic in Q 
and quadratic order in gradients [371 EH EHj ■ There are eight such terms and it is unrealistic 
to consider them all. We include one chiral L31 and two achiral terms L34 and L38. The 
L34 term was chosen as it contributes equally to all three elastic constants and the L38 term 
as it gives the largest distinction between splay and bend and because it has the largest 
contribution to the energetics of an isolated double twist cylinder. (However, the choice 
is arbitrary and other coefficients may be non-zero in any given material, and may have a 
considerable effect on the phase diagram.) 

Rewriting the free energy in terms of dimensionless variables demonstrates that it depends 
only on the two dimensionless parameters [2S] 



known as the reduced temperature and the chirality respectively, together with ratios of the 
elastic constants. In line with previous work, phase diagrams will be presented in terms of 
these parameters. 

B. Numerical minimisation of the free energy 

The free energy ([T]) is minimised by relaxing the Q-tensor according to a Ginzburg-Landau 
equation 



r := 



9(3-7) 
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This equation can be solved using many different numerical approaches; the results presented 
here were obtained using a lattice Boltzmann algorithm l¥T] . 

To study the different blue phases it is necessary to implement appropriate initial condi- 
tions for the simulation. The Q-tensor is initialised using analytic expressions appropriate 
to the high chirality limit {k oo), which act to define the symmetry of the chosen phase. 
For blue phase I [3 [23l |25] we use 



Qxx ~ ~ sm{ky/ a/2) cos{kx / V2) — sm{kx / V^) cos{kz/V2) 
+ 2 sm{kz/V2) cos{ky/V2), 

(4) 

Qxy ~ —V2sm{kx/V2) sin{kz/ V2) — V2 cos{ky / \/2) cos{kz / V2) 



+ sm{kx / V2) cos{ky/ V2), 

where k = 2\^iT/a, with a the lattice constant, and the other components are obtained by 
cyclic permutation. Similarly, for blue phase II the Q-tensor is initialised as [71 |23l EH] 

Q^^ ~ cos{kz) - cos(%), 

(5) 

Qxy ~ sm{kz), 

where k = 2tt /a and the other components are again obtained by cyclic permutation. 

Under numerical evolution using Eq. ([s]) the system relaxes to the structure of the same 
symmetry that locally minimises the free energy. We are therefore able to obtain, for 
any value of the parameters, local minima of the free energy corresponding to each of the 
cholesteric and blue phases. The global free energy minimum was taken to be the smallest 
of these calculated local minima. 

To achieve a full minimisation of the free energy it is necessary to set the correct unit 
cell size in the simulation. This is not known a priori, but rather depends on the magnitude 
of the order parameter, a quantity that is only determined by the numerical minimisation 
itself. Therefore we must introduce a means of determining, and setting, the unit cell size 
as the Q-tensor evolves during the simulation. It is possible to account for a change in unit 
cell size by rescaling the gradient contributions to the free energy and molecular field. This 
is accomplished in practice by changing the elastic constants to 

T _ T init 2 
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where a superscript 'init' denotes the initial value of a simulation parameter and r is the 
appropriate rescaling factor, which is identical to the 'redshift' described in [231 [21] ■ 

To calculate the optimal size of the unit cell we note that, since the free energy is quadratic 
in gradients, it may be written formally in k-space as 

f = ak'^ + bk + c, (7) 

where the coefficients a, b and c depend on the Q-tensor, but not on k. The optimum 
wavevector is given hj k = —b/2a, and since the coefficients a and b are determined by the 
simulation it is straightforward to use these values to determine the exact value for the size 
of the unit cell at every timestep, thereby obtaining a full minimisation of the free energy. 

III. THE EQUILIBRIUM PHASE DIAGRAM 

A. Revisiting the analytic calculations 

The phase diagram for chiral liquid crystals obtained for a selection of parameter values 
using the Landau-de Gennes free energy, Eq. ([T]), is shown in Fig. |4j Fig. |4|^a), adapted 
from [21], shows the phase diagram calculated analytically in the high chirality limit, in 
the one elastic constant approximation, and Fig. |4](b) compares numerical results for the 
same parameters in the free energy. The differences between the analytic and numerical 
results show that including higher order harmonics in the minimisation of the free energy 
is significant, as might be expected for phases with small free energy differences. In the 
numerical (exact) minimisation is stable over a larger range of parameters and at lower 
chirality values, in better agreement with the experimental phase diagram, Fig[2j Moreover, 
the regions of stability of the and O^^ textures found in the early studies [231 El] are 
shifted to unphysically high values of the chirality [22]. This result is again consistent with 
experiment, where blue phase structures of these symmetries have not been observed. 

B. Varying the elastic constants 

We now compare the phase diagrams obtained as the ratios of the elastic constants are 
varied. To investigate the effect of the bend elastic constant we chose parameter values 
L21 = L22 = -Z^34 = 0.02, L38 = which corresponds to a ratio of splay to bend of about 0.5, 
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while splay and twist remain degenerate. The resulting phase diagram is shown in Fig. |4[c). 
Comparing to the case of equal elastic constants the stability of blue phase I is seen to 
decrease quite significantly relative to the cholesteric phase while at the same time there is 
a small increase in stability over blue phase II. There is only a minor shift in the cholesteric- 
blue phase I phase boundary at the transition temperature, however, as the temperature 
decreases the shift becomes larger. 

The value of the twist elastic constant is controlled by the Landau-de Gennes parameter 
L22- In most liquid crystals the twist elastic constant is smaller than either splay or bend. 
In order to match this, we constructed the phase diagram for parameter values L21 = 
0.02, L22 = 0.04, L34 = L38 = 0, which is shown in Fig. ^d). This choice of parameters 
resulted in a ratio of splay to twist of about 1.5, while splay and bend remained degenerate. 
Again we observe that the stability of blue phase I is reduced relative to the cholesteric 
phase by an amount similar to that seen by varying the bend elastic constant. 

Finally, we consider the effect of the chiral cubic invariant on the blue phases. We chose 
parameter values of L21 = L22 = -^34 = 0.02, L^g = 0.04, L31 = —0.18, which gives a ratio of 
bend to splay of about 1.6. Moreover, for these coefficients the cholesteric undergoes helical 
sense inversion |l2l US] at a reduced temperature of about r ^ —2. The phase diagram 
is shown in Fig. |4]^e). What is remarkable is the dramatic increase in stability of blue 
phase I relative to the cholesteric phase. The region of stability has been increased down to 
chiralities as low as k = 0.07 and at such low chiralities the phase boundary is essentially 
independent of k, for all r. In addition, there is a very small region of stability for blue phase 
II located close to the isotropic transition (Fig. |4]^f)). Since blue phase I is now stable over 
a much larger temperature range it displays a significant variation in unit cell size as the 
temperature is lowered, with the lattice parameter more than doubling between r = 1 and 
r = -5. 

C. Electrostriction 

The response of the cubic blue phases to the application of an external field has been a 
topic of interest since the mid 1980s and continues to be so because of the importance to po- 
tential blue phase based devices [HI [221 HI] . The principal features include electrostriction, a 
continuous distortion of the shape and size of the unit cell with increasing field, and a series of 
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field induced textural transitions. The electrostriction involves a shift of the back-scattered 
Bragg peak of 5-10% and is quadratic in the field strength [221 EDI EI]- The direction of the 
shift changes sign with the sign of the dielectric anisotropy, but blue phase I also displays 
an unusual response referred to as anomalous electrostriction, where an expansion along the 
field direction is seen when the field is applied parallel to the [Oil] direction, but a con- 
traction for fields parallel to [001]. At larger field strengths new blue phases appear. Three 
distinct field induced textures have been identified, possessing tetragonal, screw hexagonal 
and two-dimensional hexagonal symmetry with increasing field strength [321 E31 EH ES] . 

Much of the electric field behaviour has been understood theoretically via extensions of 
the Landau-de Gennes theory, including the qualitative features of the electrostriction and 
the field induced textural transitions PSJ HEJ SZl HE], HHl [50]. However, the approximations 
inherent in the analytic calculations limited the quantitative comparison that was possible 
and a number of features, including the anomalous electrostriction of blue phase I and the 
tetragonal field induced texture blue phase X, could not be accounted for [IHl 130] . 

The main additional difficulty in minimising the free energy of the blue phases in the 
presence of an electric field [28] is in accounting for electrostriction, as there is a change 
not only in the size, but also in the shape, of the unit cell as the field is applied. This 
distortion can be accounted for in two steps: first the shape of the simulation unit cell is 
fixed and the corresponding size, which minimises the free energy, is determined. This is 
then repeated for a set of varying shapes of the unit cell, e.g. from cubic to tetragonal. 
Hence the free energy is determined for several values of the ratio L^jL^ parameterising the 
cubic-tetragonal distortion and fitted to a quadratic. The actual distortion is then given by 
the miminum of this fit. 

The field dependence of the blue phase lattice parameters is shown in Fig. [5] for an 
applied electric field along the [001] direction. Note, in particular, that the unit cell expands 
along the field direction in blue phase II, but contracts in blue phase I, in agreement with 
experiment [B]. When the field is instead applied along to the [Oil] direction both blue 
phases undergo an expansion parallel to the field, a precursor to the transition to blue phase 
X that is observed at larger field strengths. It is very pleasing that a numerical approach 
can predict anomalous electrostriction in blue phase I as the effect is lost in the truncations 
needed in analytic calculations [50]. Mapping from physical to simulation units gives a 
magnitude of the elctrostriction in the range 10^^ — 10^^/im^ V'"^ for both blue phases. 
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again in good agreement with experiments |44] . 



IV. HYDRODYNAMIC EQUATIONS 

The Ginzburg-Landau equation, Eq. ([s]), describes the relaxation of a hquid crystal to 
the minimum free energy, but does not describe physical dynamics in situations where flow 
is important. Several authors have recently used numerical approaches to solve the full 
hydrodynamic equations of motion for liquid crystals in the nematic and cholesteric phases. 
For example, it has been possible to simulate defect hydrodynamics [HH |52l [53l [M], phase 
ordering j55l |56l [57] , the kinetics of transitions between different liquid crystal phases [281 EH] ) 
cholesteric rheology [591 EOl ED ESI |63l [G^ and the effect of flow on device switching [651 
EHl EZl EHl EHl EOl ED [Z21 E3] ■ Because of their disclination structure the kinetics and the 
rheology of the blue phases is an exciting, but demanding, numerical problem which requires 
intensive numerical resources. However, such simulations are rapidly becoming feasible: we 
summarise the equations of motion and some of the results obtained so far, and then discuss 
possible directions for future work. 

The hydrodynamic equations of motion of liquid crystals are complex, both because of the 
anisotropy of the molecules and because of the coupling between the order parameter field 
and the flow fleld. In general, flow leads to a rotation of the local orientation, which in turn 
influences the flow. Similarly, if a disturbance is initiated in the director, its reorientation 
is generally accompanied by fluid motion, an effect sometimes referred to as backflow. If a 
system is close to a phase transition, or contains disclinations, variations in the magnitude 
of the order parameter can be significant and therefore a hydrodynamic description based 
on the Q-tensor is needed [H UnilZZl EH] • 

The order parameter evolves towards the minimum of the free energy, but with a convec- 
tive time derivative to account for the advection with the fluid 



The term in brackets on the right hand side is called the molecular fleld, H, and F is a 
collective rotational diffusion constant. The material derivative for rod-like molecules is 




(8) 
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given by [73] 

AQ = (5i + u ■ V) Q - (^D + n) (Q + |I) - (Q + |I) (^D - n) 
+ 2^ (Q + |I) tr (QW) , 

where D = (W + W"^)/2 and 17 = (W — W^)/2 are the symmetric and antisymmetric 
parts, respectively, of the velocity gradient tensor Wajs = V pUa- The constant ^ depends on 
the molecular details of a given liquid crystal. 

The fluid velocity field is taken to obey the continuity equation and a Navier-Stokes 
equation with a stress tensor generalised to describe liquid crystal hydrodynamics [71] 

dtQ + v ■ (f?u) = , (10) 

Q{dtU + u ■ Vu) = - Vj9 + V ■ <T , (11) 

— ^Hay{Q^p + — ^{Qa-y + \^a'y^H^i3 (12) 

5F 



Eqs. ([8]) , ( 10 ) and ( 11 ) can be solved either using a lattice Boltzmann approach [lOl [79] , or 
finite difference schemes [531 ESI |80] or, giving somewhat improved stability while retaining 
the advantages of lattice Boltzmann, a hybrid scheme where a lattice Boltzmann solution of 
the flow equation is coupled to a finite difference solution of the order parameter field 



V. BLUE PHASE RHEOLOGY 



A. Response of the disclination lattice to a Poisueille flow 

To investigate the rheological response of the blue phases we placed a unit cell between 
fixed plates and imposed constant force on the fluid, together with no-slip boundary con- 
ditions on the velocity field at the plates. In a Newtonian fluid this geometry leads to a 
quadratic, Poiseuille flow profile. The choice of boundary conditions for the director field was 
to assume that the disclinations are fixed at the boundaries. Fig. |6] compares the apparent 
viscosity (obtained comparing to the Poiseuille velocity) in blue phases I and II, in a phase 
comprising a square array of disclinations with two-dimensional cross section correspond- 
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ing to Fig. [T} and in the isotropic phase. The corresponding dischnation configurations, 
comparing zero and a finite velocity field, are shown in Fig. [7} 

For small forcing the blue phase viscosities increase by a factor ~ 4 over that of the 
isotropic fiuid. This is because the dischnation network acts to oppose the fiow and dissipate 
energy. The blue phases reach a stationary state in which the dischnation network is bent 
and twisted by the fiow. The viscosities of blue phases I and II are approximately constant 
over a range of forcing, but the square lattice structure shows shear thickening. This is 
because each defect fine, of topological strength 1, opens to a dischnation ring, comprising 
disclinations of strength 1/2. The ring then twists and bends on itself as the fiow increases, 
as shown in Fig. [7j 

As the forcing increases there is significant shear thinning in all three blue phase struc- 
tures. This occurs because the dischnation network is destroyed by the fiow, and the viscosity 
drops to that of an isotropic liquid crystal. 

These results indicate that blue phase rheology is extremely rich, and worthy of further 
study, both experimentally and numerically. There is a need for careful experiments, with 
good control over boundary conditions. On the numerical side, work is in progress to move 
towards larger numbers of unit cells in the simulations and to assess the effects of different 
flow geometries and boundary conditions [82] . 

B. Blue phase X 

The ability to solve the dynamical equations of motion means that it is possible to follow 
the way in which disclinations rearrange during continuous textural transitions. Moreover, 
the simulations provide a way to identify possible candidate structures for blue phases which 
have not yet been identifled. 

An example of this approach is simulations in which blue phase I is placed under an inter- 
mediate electric fleld applied parallel to the [Oil] direction. This is the set-up corresponding 
to the blue phase I-blue phase X transition observed experimentally in the 1980s [331 El] ■ 
The electrostriction distorts the shape of the unit cell until it becomes tetragonal, at which 
point there is a transition from the blue phase I texture to a new texture known as blue 
phase X. Starting from this geometry, and an unperturbed blue phase I, the evolution of 
the dischnation network is simulated numerically. Initially, the disclinations in the network 
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twist, they then merge to form a transiently branched structure, which finally reorganises 
into a new defect network, not stable at zero field as shown in Fig. |8j This is a candidate 
structure for blue phase X, as it is a new network, found via a continuous reorganisation 
starting from blue phase I, and only stable in a field. The results predict that in the can- 
didate blue phase X (i) the disclinations perpendicular to the field are largely unaffected, 
(ii) the network conforms to the space group identified in Reference [31], and (iii) the 
double twist cylinders deform but do not break during the transition; this observation may 
be verified by experiments along the lines of those in References [5^ IM] . 

C. Flexoelectric blue phases 

In 1969 Meyer introduced the concept of fiexoelectric coupling to an external electric 
field and showed that this could lead to a one- dimensional splay-bend distortion of the 
nematic director field [85]. These results have recently been extended to show that, near the 
isotropic-nematic transition and with sufficiently strong coupling, two-dimensional splay- 
bend structures with hexagonal symmetry can be stable |86j . 

Fig. [9] shows the evolution of blue phase I when fiexoelectric coupling to an applied field 
is increased quasistatically. As in the dielectric case, the electric field induces a twist in the 
disclination lines allowing them to transiently merge, thus fascilitating textural transitions. 
For the largest applied field strengths a transition is indeed observed to the two-dimensional 
hexagonal flexoelectric blue phase, as expected. However, at intermediate values, two further 
transitions precede it, each yielding distinct textures that are stable over a small range of 
field strengths. The first is to a centred tetragonal texture with space group J4i22. This 
has the same disclination network as the blue phase X structure found under dielectric 
switching and described above. At slightly higher field strengths there is a second transition 
to a distinct tetragonal texture with space group P4222. In this texture, the disclinations 
are all parallel to the field, occuring in pairs that wrap around each other to form a double 
helix, and with the axes of the double helices themselves then arranged on a square lattice. 
Finally, as the fiexoelectric coupling is inceased, the two members of a given double helix 
transiently merge and re-separate, allowing them to both straighten out and to adopt an 
hexagonal configuration in the plane perpendicular to the field. 
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VI. DISCUSSION 



The Landau-de Gennes expansion of the free energy of hquid crystals has proved a vital 
tool in understanding their thermodynamics. The theory has proven successful even when 
very subtle energy-entropy balances, such as those that stabilise the blue phases, come 
into play. Beautiful early calculations, using approximate theories based on the Landau-de 
Gcnncs formalism that are analytically tractable, helped understand many of the features 
of the blue phases. Here we have shown that it is possible to make further progress by 
exploiting modern computational resources to minimise the free energy exactly. 

We have shown that significant quantitative differences in the phase diagram arise from 
retaining cubic order terms in the free energy expansion. In particular, choosing expansion 
coefficients appropriate to systems where the cholesteric undergoes helical sense inversion 
gives rise to a very significant increase in the stability of blue phase 1. It would be interesting 
to see if this link can be established experimentally. Moreover, by minimising the free energy 
with respect to both the order parameter field and the size and shape of the unit cell it is 
possible to predict the electrostriction of the blue phases, obtaining good qualitative and 
quantitative agreement with experiments. 

The major drawback of the Landau-de Gennes expansion is in the difficulty of knowing 
values of the expansion coefficients for any particular compound. Inevitably this programme 
becomes increasingly difficult as additional terms arc added. Therefore the value of the 
theory is primarily in suggesting trends, and in identifying regions of parameter space where 
novel behaviour might be observed. 

Although the major features of the equilibrium behaviour of the blue phases are un- 
derstood, much less is known about their hydrodynamics. The current simulations are on 
tiny samples and it is necessary to guess suitable boundary conditions on the director field. 
However, the rcponsc of the blue phases and their disclination networks to an imposed 
flow is fascinating, and a programme of developing the numerical approach in tandem with 
experiment is likely to uncover novel physics. 

Solving the equations of motion of the blue phases also makes it possible to investigate 
their dynamics under changes in an applied electric field. Understanding this is an important 
step to assisting device design. Moreover it is possible to follow the kinetics of transitions 
between the blue phases, as an external parameter such as the electric field is varied. This 
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allows prediction of possible sequences of phases that will, in general, depend on the values 
of the free energy, whether there is a convenient path through phase space allowing any 
given phase to be accessed, and the speed with which the perturbation is applied. For the 
examples described here the structural evolution depended primarily on relaxation kinetics, 
with hydrodynamics playing a minor role. However, this may not always be the case. 

Finally we return to blue phases that have been stabilised over a larger temperature range. 
One approach has been to use bimesogenic molecules, and there remain interesting questions 
as to whether the resulting blue phase is stable or metastable and in identifying the physical 
mechanisms behind the increased temperature range. The polymer stabilised blue phases 
can be interpreted in terms of the polymers pinning the disclinations. As numerical work 
on ever more complex fluids develops, there are many exciting questions about blue phases 
in mixtures of liquid crystals and polymers, colloids and nanoparticles, that will become 
increasingly accessible. 
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FIG. 1: Director field in a hypothetical, two-dimensional, blue phase showing how local regions of 
double twist can be pieced together with a square array of topological defects. (Blue phases with 
two-dimensional symmetry and translational invariance in the third dimension have been observed 
experimentally in an electric field [32|, but these have hexagonal symmetry.) 
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FIG. 2: Experimental phase diagrams of the cholesteric blue phases for two different chiral com- 
pounds, reproduced from Reference [4 . Three distinct blue phases are found, in the order blue 
phase I, blue phase II, blue phase III upon increasing the amount of chiral dopant. 
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tb) 

FIG. 3: Network of disclinations and arrangement of double twist cylinders in blue phase I (a) and 
blue phase II (b). The disclinations are shown on the left in blue and the double twist cylinders 
on the right in grey. In all figures 2^ unit cells have been shown to more clearly illustrate the 
structure. 
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FIG. 4: Phase diagram of cholesteric liquid crystals within the Landau-de Gennes theory, r, the 
reduced temperature, and k, the chirality, are defined by Eq. Top row: phase diagrams in 
the one elastic constant approximation determined (a) analytically using truncated Fourier series, 
and (b) numerically. Middle row: numerical phase diagrams with unequal elastic constants: (c) 
Ki = K2 = O.S-fTs, (d) Ki = Ks = 1.5K2. Bottom row: (e) the numerical phase diagram obtained 
with the chiral invariant, L31, added to the free energy. The magnitude of this term was chosen 
so as to produce helical sense inversion in the cholesteric phase at a temperature not far below 
the isotropic transition temperature, (f ) an enlargement of the region near the isotropic transition 
temperature. Note the reversal in the order of appearance of blue phase I and blue phase II as a 
function of chirality. (adapted from p6] ) 
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FIG. 5: Distortion of the unit cell of blue phase I (top) and blue phase II (bottom) in an electric field 
applied along [001]. The lattice parameters parallel (o, solid line) and perpendicular (□, dashed 
line) to the field are shown relative to their zero field value. Note that blue phase I contracts 
parallel to the field, whereas blue phase II expands parallel to the field, (adapted from [28] ) 
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FIG. 6: Variation of the apparent viscosity of blue phase I, blue phase II, and a square array of 
double twist cylinders as a function of the applied forcing. The viscosity of the isotropic liquid 
crystal is also shown. Initially the blue phases show an enhanced viscosity as the disclination 
network opposes the flow. At higher forcing the disclination network breaks up and the viscosities 
of the blue phases tend to that of the isotropic liquid crystal, (adapted from |27| ) 
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FIG. 7: Structure of the blue phases under Poiseuihe flow. Rows A, B and C correspond to a square 
array of double twist cylinders, blue phase I, and blue phase II respectively. The first column shows 
the disclination network at zero forcing and the second a steady state of the network under flow. 
Note how, in row A, the strength 1 disclination opens to a strength | ring, (adapted from [27J ) 
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FIG. 8: Simulation of the blue phase I-blue phase X transition. Top: blue phase I before the 
field is applied. Middle: at the transition from blue phase I to blue phase X. Bottom: the final 
configuration of disclination lines in blue phase X. In all cases the electric field is applied parallel 
to the [Oil] direction, which is vertical in the figure, (adapted from |^5] ) 
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FIG. 9: Textural transitions in blue phase I induced by flexoelectricity. The electric field is applied 
parallel to the [Oil] direction (vertical in the figure) and produces a series of transitions first to 
two distinct textures possessing tetragonal symmetry, with space groups /4i22 and P4222, and 
finally to a two dimensional hexagonal texture (viewed here along the field direction) at larger field 
strengths. 
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